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■ Abstract 
[JL ■ 

'^i ■ We consider spatial population dynamics given by Markov birth-and- 

■ death process with constant mortality and birth influenced by establish- 
, ment or fecundity mechanisms. The independent and density dependent 

dispersion of spreading are studied. On the base of general methods of [14j . 
we construct the state evolution of considered microscopic ecological sys- 
tems. We analyze mesoscopic limit for stochastic dynamics under consid- 
^Sj ■ eration. The corresponding Vlasov-type non-linear kinetic equations are 

^ ' derived and studied. 
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1 Introduction 

Complex systems theory is a quickly growing interdisciplinary area with a very 
broad spectrum of motivations and applications. One may characterize complex 
' systems by such properties as diversity and individuality of components, local- 

ized interactions among components, and the outcomes of interactions used for 
replication or enhancement of components. In the study of these systems, proper 
language and techniques are delivered by the interacting particle models which 
form a rich and powerful direction in modern stochastic and infinite dimensional 
analysis. Interacting particle systems are widely used as models in condensed 
matter physics, chemical kinetics, population biology, ecology, sociology, and 
economics. 
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Mathematical realizations of such models may be considered as a dynamics 
of points in proper state spaces. In some applications the possible locations for 
the points of system are structured, e.g., if we consider dynamics on graphs, or, 
in particular, on lattices. Another class of models can be characterized by the 
free positions of points in continuum, e.g., in Euclidean space M''. As it was 
shown originally in statistical physics, many empirical effects, such as phase 
transition, are impossible in systems with finite number of points. Therefore, 
systems with infinite points can be considered as mathematical approximation 
for realistic systems with huge but finite number of elements. Among all infinite 
systems we will study locally finite ones. Namely, the configuration space over 
space M.'^ consists of all locally finite subsets (configurations) of M.'^ 

r = r{R'^) := {7 C M.'^ |7a| < 00, for all A e 6b(K'')}. (1-1) 

Here 7a := 7!"! A, the symbol | • | stands for the cardinality of a set, and /i3b(R'^) 
denotes the class of all bounded Borel sets in R'^. Each configuration may be 
identified with a Radon measure on R** by the relation 7(A) — \ja\- As a 
result, r can be equipped with the vague topology and the corresponding Borel 
(T-algebra. 

Depending on application the points of system may be interpreted as molecules 
in physics, plants in ecology, animals in biology, infected people in medicine, 
companies in economics, market agents in finance, and so on. It is supposed 
that points of a system evolve in time interacting with each other. In the present 
paper we focus our attention to the dynamics with birth and death mechanisms. 

The spatial birth-and-death dynamics describe an evolution of configura- 
tions in R'*, in which points of configurations (particles, individuals, elements) 
randomly appear (born) and disappear (die) in the space. Heuristically, the 
corresponding Markov generator has the following form: 



{LF){^)=Y,d{xn\x)D-F{^)+ [ b{x,^)D+F{^)dx, 



(1.2) 



a;G7 

where for F : F ^ M, x ^ 7 

D-F{^)^F{^\x)~F{^), D+F{^) = F{^Ux)-F{^). (1.3) 

Here functions d and b describe rates of death and birth correspondingly (for 
details see, e.g., [H]). 

In the present paper we apply the results of [13] to study the question about 
the existence of the evolution corresponding to (|1.2p for a particular choice of 
the functions d and b. This question can be answered once we will be able to 
construct a semigroup associated with L in a proper functional space. This 
semigroup determines the solution to the Kolmogorov equation, which formally 
(only in the sense of action of operator) has the following form: 

= LFt, Ft Fo. 



dt 
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To show directly that L is a generator of a semigroup in some reasonable func- 
tional spaces on F seems to be difficult problem. This difficulty is hidden in 
the complex structure of non-linear infinite dimensional space T. However, in 
various applications the corresponding evolution of states (measures on con- 
figuration space) helps already to understand the behavior of the process and 
makes possible to predict the equilibrium states of our system. In fact, proper- 
ties of such an evolution itself are very important for application. The evolution 
of states is heuristically given as a solution to the dual Kolmogorov equation 
(Fokker-Planck equation) : 

— ^ L nt, nt |j^(,= fJ.0, (1.4) 

where L* is an adjoint operator to L defined on some space of measures on F, 
provided, of course, that it exists. 

Technically, we will study solutions of (|1.4p in terms of correlations functions, 
n >0 which are symmetric functions on (K'')" and related to a density of 
distribution for each n points of our system (rigorous definition will be given in 
Section 2). 

Among all birth-and-death processes we will consider only those in which 
new particles appear from existing ones. These processes correspond to the 
models of the spatial ecology. In the recent paper [T^], we studied Bolker- 
Dieckmann-Law-Pacala ecological model, which corresponds to the following 
mechanism of evolution. Each existing individual can give birth to the new 
one independently of all other individuals of the system. It may also die in- 
fluenced by the global regulation (mortality) again independently of all other 
members of the population or it dies because of the interaction with the rest 
of the population (local regulation). The latter mechanism may be described 
as a competition (e.g., for resources) between individuals in the population. 
Heuristically, the corresponding Markov generator has the form (|1.2|) with 

d{x,^) ~ 771 + >f:^ (a; — y), (1-5) 

h(x,i) = >f+^a+(a; - y), (1.6) 

Here a+,a~ are probability densities, and constants m^>i^^>t' > 0. In popu- 
lation ecology, the constant 777 is called mortality and the functions a+,a~ are 
known as dispersion and competition kernel, respectively. 

By [12], if 777 = >c^ = (free growth model) then the first correlation 
function (density of the system) grows exponentially in time. To suppress this 
growth we may consider the case 777 > >c~ — (contact model, see also [T51I20] ). 
Then for 777 > >c~^ we obtain globally bounded density (even decaying in time 
for m > ><+). Nevertheless, locally the system will show clustering. Namely, 
/Cj""* n! on a small regions for t > (see [12] for details). The main result 
of [12j may be informally stated in the following way: if the mortality m and 
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the competition kernel yt a are large enough, then the dynamics of correla- 
tion functions associated with the pre-generator (|1.2p preserves (sub-)Poissonian 
bound for correlation functions for all times, i.e., fe^"'' < C", C > 0, n > 1. 

In the present article we introduce new mechanisms of local regulation in the 
corresponding system, alternatively to (|1.5p . Namely, we set >c~ = in (jl.Sp 
and consider two different modifications of (jl.6p . The first one includes the 
influence of the whole system on the reproduction (fertility, fecundity) of each 
single individual. The second modification of (jl.6p contains a mechanism which 
shows establishment of each individual in the system. The precise descriptions 
are given in the next section. Such models have been actively studied in modern 
ecological literature, see e.g. [S] and references therein. Here, for the first time, 
we present a rigorous mathematical description for these evolutions. 

This article is organized in the following way. In Section 2, we describe the 
model rigorously providing the proper spaces for the corresponding functional 
evolutions. In Section 3 we apply general results about birth-and-death dynam- 
ics on configuration spaces obtained in |14j . Informally, the main results state 
that if mortality m is big enough and negative infiuence of establishment or 
fecundity is dominated by dispersion then the corresponding evolution exist. In 
Section 4, we study the mesoscopic description of our model in terms of Vlasov 
scaling. 

It should be noted also, that the Vlasov-type scalings for some Markov pro- 
cesses on finite configuration spaces were considered in (2HS]- Note that the 
corresponding limiting hierarchy was obtained at the heuristic level. In the 
present paper, we prove a weak convergence to the limiting hierarchy in the 
case of infinite continuous systems for bounded but non-integrable densities. 

It is worth pointing out that the necessity of a big mortality is a result 
of perturbation theory for linear operators which gives the existence of the 
corresponding dynamics for the infinite time interval. However, with the help 
of another technique considered in [7] , [10] , we are able to show the existence of 
the dynamics with any mortality but only on finite interval of time. This result 
will be presented in the forthcoming paper. 

2 Description of model 

We recall that the configuration space F is given by (|l.ip . It is equipped with 
the vague topology, i.e., the weakest topology for which all mappings F 3 7 i->- 
'^xe-y /(^) S ^ ^"^^ continuous for any continuous function / on M'' with com- 
pact support. The space F with the vague topology is a Polish space (see, 
e.g., [16] and references therein). The corresponding Borel cr-algebra S(F) will 
be the smallest cr-algebra for which all mappings F 3 7 n- |7a| G Nq := N U {0} 
are measurable for any A e Bh{^'^), see, e.g., ilj. We set J^;yl(F) for the class of 
all cylinder functions on F. Each F G J^;yi(F) is characterized by the following 
relation: ^(7) = F(7a) for some A e BbiR"^). 
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Let < ^ € L^{R''') be given even function such that 



C0 



(l -e-^^^^^dcc G (0; +00). (2.1) 



For any even < / € we denote 

EHv) := E E /(^ - 2^)' e 
Ef{x, 7) := E - 7 e r, X G 

if) := [ fix)dx. 

As it was already mentioned in the Introduction wc would like to stiidy two 
classes of the interacting particle systems (IPS), whose mechanisms of evolution 
are described by the corresponding heuristically given Markov generators: 



(LestF)(7) := m^[F(7\x)-F(7)] 

+ E / M^,yn\ y)e-^*(^'^) [F(7 U ar) - ^(7)] rfa(2.2) 



and 

(Lfeci^)(7) := m5^[F(7\ar)-F(7)] 



+ V e-^'(^''^\^) / 6o(x, y, 7 \ y) [i^(7 U x) - i^(7)] d^^-S) 

The first model shows the influence of establishment in the system and the 
second one presents fecundity. Here and in the sequel the mortality m is al- 
ways supposed to be strictly positive. One can see that the establishment rate 
g-E (a;, 7) -^iii be smaller if x will be inside or close to the dense region of the 
configuration 7. In its tiun the fecundity rate e"-^'''^^''''^^' would be also smaller 
if y is situated in the dense area of 7. The non-negative measurable rate 60 
represents the dispersion of the model. Let < a+,6+ G L^(]R'') be given even 
functions, and (a+) = 1. We consider two types of the dispersion: 

• density independent dispersion 

6o(a;,y,7\y) = x+a+{x - y), 

• density dependent dispersion 

bo{x,y,-f\y) = a+{x-y)(H++ ^ b+{y-y')]. 

y'^i\y 
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As it was mentioned above, we will study evolution of our model in terms 
of its correlation functions. Below we introduce some basic notions needed to 
describe the corresponding evolution. 

The space of n-point configurations in an arbitrary Y G B{M.'^) is defined by 



r(")(r) := I?] C r |?7| = n|, nen 



By definition we take T^°\Y) {0}. As a set, r(")(r) may be identified 
with the symmetrization of F" = {(xi, . . . ,x„) G \ Xk ^ Xi if k ^ l}. 
Hence one can introduce the corresponding Borel cr-algebra, which we denote 
by S(r(")(r)). The space of finite configurations in an arbitrary Y € ^(M'') is 
defined by 

ro(y) □ r(")(y). 

This space is equipped with the topology of the disjoint union. On ro(y') we 
consider the corresponding Borel cr-algebra denoted by B(ro{Y)) . In the case of 
F = M'^ we wiU omit Y in the notation. Namely, Tq := ro(R''), F'") r(")(]R''). 

The restriction of the Lebesgue product measure (dx)" to (r("\ ;B(r("^)) we 
denote by m^"). We set m^°^ := S^^y The Lebesgue-Poisson measure A on Fq 
is defined by 

n=0 

For any A G Bhi'R.'^) the restriction of A to F(A) := Fo(A) wih be also denoted 
by A. The space (F,S(F)) can be obtained as the projective limit of the family 
of spaces {(F(A),i3(F(A)))}^gg^^jgj^, see, e.g., [1]. The Poisson measure tt on 

(F, B(r)) is given as the projective limit of the family of measures {7i''^}Aeeb(R'')> 
where tt^ := e~™^'^''A is the probability measure on (F(A), S(r(A))) and to(A) 
is the Lebesgue measure of A G Sb(M''); see, e.g., [T]. 

A set M G S(Fo) is called bounded if there exists A G i3b(M'*) and N G 
N such that M C U^=o ^'■"H^)- The set of bounded measurable functions 
with bounded support we denote by i?bs(Fo): i-e-i G G i3bs(Fo) if G \ro\M= 
for some bounded M G B{To). Any B(Fo)-measurable function G on Fq, in 
fact, is defined by a sequence of functions {C*^"-* j^^j^ where G*^"-' is a S(F'"^)- 
measurable function on F^"). As usual, functions on F are called observables 
and functions on Fq are called quasi-observables. 

There exists a mapping from _Bbs(Fo) into 7\;yi(F), which plays the key role 
in our further considerations. It has the following form 

ifG(7):=^G(r;), 7 G F, (2.4) 

where G G Bbs(Fo), see, e.g., [l5|2H22j . The summation in ()2.4p is taken over all 
finite subconfigurations 77 G Fq of the (infinite) configuration 7 G F; we denote 
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this by the symbol, 77 d --y. The mapping K is linear, positivity preserving, and 
invertible, with 

K-'F{i^):=Y,{-lt\i\F{0, rjeV^. (2.5) 

Note that if function F has special form 

F(7)-^i?(a;,7\a:), 

where H{x^ ■) is defined point-wisely at least on Fq, then, by direct computation, 
{K-^F){i^)^Y.i^'^H{x,-)){ii\x), r/GTo. (2.6) 

We set also 

{K,G){ii) := {KG){i^), rye To. 

Let us define, for any ;B(]R.'')-measurable function /, the so-called coherent 
state 

ex{f,v)-=Y{f{^), ^?ero\{0}, eA(/,0):-l. 

Then 

(XoeA(/))(ry)=eA(/ + l,r7), e Tq (2.7) 
and for any / G L'^{W^, dx) 

eA(/,?7)rfA(r,) -expf / f{x)dx]. (2.8) 

A measure ^ S 7Mfj^(r) is called locally absolutely continuous with respect 
to the Poisson measure tt if for any A G ,Bb(IR.'*) the projection of /i onto r(A) is 
absolutely continuous with respect to the projection of tt onto r(A). By [15], in 
this case, there exists a correlation functional fc^ : Fq — > such that for any 
G G i?bs(ro) the following equality holds 

(XG)(7)d/i(7) = / G(77)fc^(77)dA(r;). 

(n) . . in) 

The restrictions fc/j of this functional on Fg , n G No are called correlation 
functions of the measure ^. Note that k'^^'' = fc^(0) = 1. 

We recall now without a proof the partial case of the well-known technical 
lemma (see e.g. llQp which plays very important role in our calculations. 

Lemma 2.1. For any measurable function iJ : Fq x Fq x Fq — > M 

/ 5]i?(e,'7\e,'7)rfA(r7)= / / H(e,ry,77UOdA(OdA(r,) (2.9) 

if both sides of the equality make sense. 
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For arbitrary and fixed C > 1 we consider the functional Banach space 

Cc.= L\To,C\^\\{di^)). (2.10) 

In the sequel, symbol ||-||p stands for the norm of the space (|2.10p . 
Let d\c := C' ldA, then the dual space 

{Cc)' = {L\To,d\c))' = L^{Ta,d\c). 

The space {Cc)' is isometrically isomorphic to the Banach space 

/Cc:={fc:ro^M fc C"! ! G L°°(ro, A)} 

with the norm 

:= ||C-l-lfc(.)|U~(r„,A) 
where the isomorphism is provided by the isometry Rc 

[Cc)' 3k^ Rck := fcCll £ ICc- 

In fact, one may consider the duality between the Banach spaces Cc and 
Kc given by the following expression 

((G, fc)) / G-kd\ GeCckelCc (2.11) 
"'ro 

with |((G, < \\G\\c ■ WkWjcc- It is clear that k 6 K-c implies 
\kir])\ < \\k\\Kc C'''' for A-a.a. 77 e Tq. 

In the paper [17] , it was proposed the analytic approach for the construction 
of non-equilibrium dynamics on F, which uses deeply the harmonic analysis on 
configuration spaces. By this approach the dynamics of correlation functions 
corresponding to (|1.4|) is given by the evolutional equation 

dkt 

= L'^ku h ko, (2.12) 

where L'^ is a dual operator to the X-image of L defined by the expression 

L R-^LK 

with respect to the duality ()2.1ip . Hence, — L* . In order to construct the 
evolution of correlation functions we are going to follow such a scheme: we show 
that L is a generator of a Go-semigroup in the certain Banach space and after 
consider the dual semigroup which solves the Cauchy problem (j2.12p . 
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3 Functional evolutions 



Let 

V:^{GeCc\ \ -\G{-)^Cc]. 

Note that -Bbs(ro) C T). In particular, 2? is a dense set in Cc- 

In [Tl], we have found sufficient conditions for operator {L,T>) to be a gen- 
erator of a semigroup in Cc- In the case of Markov generators (|2.2p or (|2.3p . 
this result may be formulated in the following way. 

Lemma 3.1 (Theorem 3.2 of [14]). Suppose there exists < a < ^ such that 



(3.1) 



(3.2) 



J2 I \{K^'h{x,{^\x)U-)){r^)\G\^\d\{i^)<am\il 

where h{x, rj) is equal either 

e-^*(^''')^6o(x,2/,ry\y) or ^ e-^*(^'''\^)&o(a:, y, ^7 \ y)- 

Then {L,T>) is the generator of a holomorphic semigroup in Ce- 
lt is worth noting that if p.ip is valid, then for any G E V 

{LG) (r,) = -m|ry|G(ry) + ^ / G(e U x) {K^'h{x, - U 0) (?/ \ Odx 

Theorem 3.2. Let < a+,6+,(/) G L^iW^) be even functions such that ([21]) 
holds and (a+) ~ 1, B := (6^) > 0. Suppose, additionally, that there exist 
constants Ai, A2 > such that 

0<a+{x) < Ai(l}{x), xeW^, (3.3) 

a+{x-y)b+{y-y') < A20(x - 2/)0(a; - y'), x,2/,y'eM^ (3.4) 

A1X+ AA2 AiB , ^2(0) ^ ^ / N 

+ ^ + + K+ + + GB < -e-"'^^ - 3.5 

eC e-^G e e 2 

Then (j3.ip /loWs a77.rf (Lost = K^^L^stK^T)) is the generator of a holomorphic 
semigroup C/ost(i) "in ^c- 

Remark 3.3. In the density independent case, 6+ = 0, the assumption p.4p 
holds with A2 — 0. Moreover, since B — 0, the condition p.5p will have the 
following form 

^ + < -e 

Before proof of Theorem 13.21 we give an example of a^, which satisfy 
(13.411 in the Lemma below. 
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Lemma 3.4. Suppose that there exist constants Ei^ E2 > and S > d such that 

E E 
a+{x) < b+{x) < ^ < E2<j)ix), X e R'^. 



Then holds with A2 ^ E^. 

Proof of Lemma \3.4\ Using obvious inequality 

l + \x-y'\<l + \x-y\ + \y-y'\<{l + \x- y\) {l + \y- y'\) 
we obtain that 

a+ix - y)b+ {y - y') < < Ej<f>{x - y)<P{x - y'), 

{l + \x-y\) {l + \x-y'\) 

that proves the statement. □ 

Proof of Theorem \3."A Let us set 



6e.t(x,7) = e-^*(--^)^a+(a;-2/)U++ ^ h+{y-y')\. (3.6) 
yei ^ v'ei\y 

To check (|3.ip . we wiU try to estimate the integral 

|(Xo-i5e,t(x,eU-)) W|Cl"ldA(r7), ^ e To 

To 

uniformly in x G M'* and ^ G Fg. 
In view of one has 

6est(a;,CUr;) = e-^*(-'«)e-^*(-''')^a+(x-y)('x++ ^Hy-y'] 

ye? y'ei\y 
^^-eHx..O^-eH-,v) ^ ^ a+(a; - y)h+{y - y') 
y'eri yei 

+^-eU^.S)^-eH-,v) J2a+{x~ y') L+ + ^ b+{y - y') 
y'ev ^ ye? 

+e-^*(-'«)e-^*(^''')^a+(x-y) ^ b+iy~y'). 

yen y'&Av 
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Using (P3|) - (PT7)) . wc obtain 



(3.7) 



e\ e 



-0 X- 



y'er] ye(, 



J2 ex (e-^(— ) - 1, \ y') a+{x - yOe-^'""^') 



v'&v 



(e-^(^-)-i,77\{y,y'}). 

Next, let K = e"*^ then, by (g^, 

< K6est(x,0 + Ce-^'^"'«^ / f y2a+{x-y)b+{y-y')e 

Jr„ Jw 

y-ex (|e-"*("-) - l| Cl^ldy'dA (77) 

(I 



-4>{x-y') 



+Ce" 



,77) a+{x-y')e-'l''^'''~y'^ 



x[>,+ + J2b+{y-y'))c\''^dy'dX{7j) 



xeA 



i+(x - y)b+{y - y')e-^(^-?')e-'^(^-^') 
l| C^'^^dy'dydX (tj) 



ye? ' 



+KC2e-^**(^'«) / / a+ix-y)b+iy~y')e-'''^''-y^e-'''^'-'-y">dy'dy. 



11 



By p.3p . one has 

e-^'(-^i)y a+{x ^ y) < A^e-'''^^'^^E^{x,0 < 

^-^ e 

where we used the elementary inequality xe~^ < e~^, x >0. Next, by p.4p . we 
may estimate 

f a+ix - y')e~'^'^''-y'^b+{y - y')dy' 

Moreover, yield 

y<^i y'&^\y 

e e e 

since x^e^"^ < 4e^^, x > 0. 
Therefore, we have 

\{K^%st {x,^L>-)) (77)|Cl''ldA {tj) 

< k{— + ^] + kCB— + kC>c+ + kC^-^ + kC^B =: D. 

\ e e'' J e e 

To obtain (I3.ip . it is enough to suppose that D < am, where % < \- Hence, we 
need that m > ^ only, that is p.5p . The theorem is proved. □ 

Theorem 3.5. Let < a+,6+,(/) e ^^(M'') &e ewen functions such that p.ip 
fto/ds anc? (a+) = 1, i? = > 0. Suppose, additionally, that there exists 

constants Ai, A2 > such that for a. a. x, y, y' G 

Q<a+{x) < Ai0(x)e-'^(^\ (3.8) 
h+{x) < A2C^{x), (3.9) 

.^ + ^ + CB+(^+B)^ + l^C<^e-^^. (3.10) 
e V C / e 2 

T/iera p.ip /lo/ds and (Lfcc = K^^LfccK,T>) is the generator of a holoraorphic 
semigroup Ufedt) in Cc- 
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Remark 3.6. In the density independent case, A2 = B = 0, and p.lOp may be 
rewritten in the form; 

Proof. Set 

6fec(x,7) = ^e-^*(''''^\^)6o(x,y,7\2/). 

Then, one has 

y6»? v'&v\y 

y'Gi; yi^i 

and, using (I2.5p - (l2.7p . we obtain 

(X„-i6fcc(a;,CU-))(r/) (3.11) 
= E e-''''^'«^eA (e-^(^-) - 1, ,7 \ y) - y) (>^+ + E ^^(^^ " 



+ Y^e~^*^y'^^a+{x-y) ^ - y')e"*(^"^')eA (e-^'^-' - 1, 77 \ y \ y') 

ye? ^ y'6«\y ^ 

^ ex (e-*(?^-) - 1, 77 \ y') e-'^(^-^') ^ e-^*(^'«\^)a+(:E - y)6+(y - y'). 
y'ev yei 



Therefore, for k = e'^*'^ we have, by (j2.8p . 

/ |(i<ro-i6fec(x,eU-))(?])|Cl''ldA(r;) 
< kC [ e-^^^y-^^a+{x^y)(>c+ + y b+{y-y')]dy 

+kC^ [ I e-^*^y^^^a+{x-y)h+{y^y')e-'l'iy-y')dydy' 
ye? ^ y'e?\y 
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+kC / e-'f'iy-y') V e~^*(«^«\^)a+(x - y)b+{y - y')dy' 

< kC>c++kc[ a+{x-y)e-^''^y'^'>y2b+iy-y')dy 
+kC2B + {kk' + kCB) ^ e-^*(f'«\2/)a+(a; _ y) 

a'e?\y 



By 



AcC / a+ix - y)e-^'(^'«) V - y')dy < "^"^^ 



Now we will verify the sufScient condition (I3.ip . By p.8p . one has 

and, by p.8p and p.9p we have 

= ^ e*(^-«)a+(x-2/) ^ b+{y-y') 

ye? 2:6{\y y'e(C\a)\i: 

y&i xei\y y'ei\y 

Hence, by (PIUI) . we obtain ([311]). □ 

Let (L', Doni(i7')) be an operator in (Cc)' which is dual to the closed op- 
erator [L^V). Here and below L means either Lost or Lfoc- We consider also 
its image on ICc under the isometry Rc, namely, let L* = RcL' Rq-i with the 
domain Dom(L*) = RcBom{L'). 

By Proposition 3.5 of [14, for any a G (0; 1) 

ICac C Dom(L*). 

Under the conditions of Theorem l3.2l or Theorem I3.5[ there exists a e (O; 
such that (13.11) holds. In the following let T{t) denotes either Ucstit) or U[cc{t). 
One can consider the adjoint semigroup T'{t) in (Cc)' and its image T*{t) in 
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fCc- By, e.g., Subsection II. 2. 6 of [2], the restriction T®{t) of the semigroup 
T*{t) onto its invariant Banach subspace Doni(L*) (here and below ah closures 
are in the norm of the space JCc) is a strongly continuous semigroup. Moreover, 
its generator L® will be part of L*, namely, 

Dom(L®) |fc e Dom(L*) L*k G Dom(L*)| 
and L*k = L®fc for any k e Dom(L®). 



Theorem 3.7 (Theorem 3.8 of [H]). For any a G (77;!) the set JCac is a 
[t) -invariant Banach subspace of ICc- 

Therefore, for a G (^; l), one can consider the restriction T®" of the semi- 
group T® onto ICac- This restriction will be strongly continuous semigroup 
with generator L®" which is restriction of L® onto ICaC (see, e.g.. Subsection 
II.2.3 of |9]). Therefore, 

Dom(L®") = |/c e L*k£lC^Y 
and L®" coincides with L* on Dom(L®"). Note that for any k G IC^c C D{L*) 

{L*k)irj)^~m\rj\kirj)\+J2 [ U iv\ ^)){Ko'bix, ■ U v\ x))iOdX{0- 

The explicit expressions can be found using p.7p or p. lip . 

Hence, we have the strong solution (in the sense of the norm in K-c) of the 
evolution equation 

?rkt = tkt (3.12) 
ot 

at least on the subspace K-ac- 

Remark 3.8. To study stationary equation L*k = corresponding to (j3.12l) on 
the set of functions k G ICaC such that fc(0) = 1, one may consider even weaker 
assumptions without denominator 2 in p.5|) or p.lO|) . However, by Proposition 
3.9 of [13], a unique solution of this equation will satisfy k{ri) — for all \t]\ ^ 0. 



4 Vlasov scaling 

To begin with, we would like to explain the idea of the Vlasov- type scaling. 
The general scheme describing this scaling for the birth-and-death dynamics 
as well as for the conservative ones may be found in [13j . This approach was 
successfully realized for the Bolker-Dieckmann-Law-Pacala model (|1.2p - (|1.6p 

in [n]- 

Let us now detail how we proceed to organize the Vlasov-type scaling. We 
will initially scale the generator L by the scaling parameter e > 0, in such a way 
that the following holds. First of all the if-image of the rescaled operator 
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has to be a generator of a semigroup on some Hc^ ■ Consider the corresponding 
dual semigroup T*{t). Let us choose an initial function of the corresponding 
Cauchy problem depending on e in such a way that k^o\v) £^l''lro(?y), e 0, 
7y e Fq with some function rg, independent of e. Secondly, the scaling L i-^ 
has to be performed to assure that the semigroup T*{t) preserves the order of 
the singularity: 

Moreover, the dynamics ro n- rt should preserve coherent states. Namely, if 
ro(ri) = e\{po,r]), then rt{r]) — ex{pt,ri) and there exists explicit (nonlinear, in 
general) differential equation for pt: 

d 

g^Ptix) = V{pt){x) 

which is called the Vlasov-type equation. 

Below we realize this approach for the case of 

(LF)(7) = mJ2D-F{^) + f b{xn)Dt F{-f)dx, 

where b — fe(a+, 6+, 0) is either birth rate with establishment (see ()2.2p ) or the 
one corresponding to the fecundity mechanism. Let us consider for any e G (0; 1] 
the following scaling 

(L,F)(7)-m^i?-F(7)+£-i / b,{x,^)D+ F{^)dx, 
xe-y -^K" 

with be = b{ea~^ , eb"^ , ecj}) . Here are given by ()1.3p . We denote by &e,est 
and 6e,fec the scaled rates for the corresponding models. We define also the 
renormalized operator (see [TT1[T3] for details) 

where {R„G){rj) = (tI''IG(77) for arbitrary cr > 0. 

Lemma 4.1. Suppose that the conditions of Theorem \3.2\ (or Theorem \3.5\) are 
satisfied with {4>) instead of in p.Sp (in p.lOp . correspondingly). Then there 

exists a £ ^0; ■^'^c/i that 

J2 [ \{K,\{x,{(\x)U-)){v)\e-\^\c\^\dX{7^)<am\^\, (4.1) 
xei-'^o 

where b^ = b^^cst (or b^ — be^fcc, correspondingly) . 
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Proof. We begin with the estabUshment case. Set 

=e-i(e-"'^(^' -l), x^R"^. 

By (13.71) . we have 

£-|''l(Xo-i6,,est(x,CU-))(77) (4.2) 
+ee-^^*(-'«) E E - - y')e-^^^"-^''eA -•),'? \ 2/') 

XGA (■0e(a;--);^\{y.y'})- 
Since e £ (0; 1] and 

\i'e{x)\ < (j){x), X e R'', 

the estimate for e^l^'l \KQ^b^^cBt ix,^LI ■) \ (77) wiU be ahuost the same as for 
li^Tp^^fe (x, C U •) I (77) in the proof of Theorem 13.21 The changes wih concern the 
term \e~'^ — 1| which will be substitute by (j)- This leads to the new constant 
(cj)) instead of in further estimates. The rest part of the proof is the same as 
for the non-scaled case. 

The same approach may be used for the case of fecundity. Indeed, 

e-l"l(i^o-i6e,fcc(x,^U-))(?7) (4.3) 

= e^'^'^^^'^^eA (My -■),V\ y) a+{x - y) L+ + ^ eb+{y - y')) 

+ ^e-^'^'(^'«)a+(x-2/) b+iy -y')e-'*^y~y')e^{My - ■),v\y\y') 

yev y'eri\y 

+sYexiMy - ■),v)e-''''^y-^'^y^a+ix - y)iK+ + ^ eb+iy-y')' 
y6« y'ei\y 

+e Y eA (My --IvX y') e-^y^y') Y e--^'^y^^\y'>a+(x - y)b+{y - y'). 
y'^ri yei 

The analogous arguments to establishment case complete the proof. □ 

Under conditions of Lemma |4. II we have the following result about the renor- 
malized semigroups in Cc and JCc- 
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Proposition 4.2 (Proposition 4.1 of [2]). Let the conditions of Lemma \4-l\ 
hold. Then for any e € (0;1], (Le^est.renj and (l/£,fec,ron, 2^) are the genera- 
tors of holomorphic semigroups C/e,est(i) and Us,iec(t) on Cc, correspondingly. 
Moreover, there exists ao G (0; i) such that for any a £ (ao; ■^) and e G (0; 1] 
there exist a strongly continuous semigroups Uf^{t) on the space fCaC with gen- 
erator L®? = L* a on the domain 



Dom(L®^) = {A: e /C„c | ^^j^rcn^ £ ^^c]. 
Here and below '(j ' means 'est ' or 'fee correspondingly. Note that, for k e JCac 
(L:,j,,,„/c)(77) - -m|77|fc(,7) (4.4) 
+ E / U ('? \ ^))^"'^' (^o"'^-.»(^' • U 77 \ x)) (OdA(0. 

By (|4.2I) . (14. 3p . there exist the fohowing point-wise limits 

lime-\^\{K^\.,,tix,^U-))iv) (4.5) 
= x+ E eA(-0(a;- •)j'7\y')a+(a;-y') 

y'eri 

wsi; y'ev\y 

and 

hme-l"! (i^o-i6,,fec (a;, •))(??) (4.6) 
= ><+EeA(-0(2/--):'7\2/)a+(x-?/) 

+ Ea+(x-y) b+iy -y')exi-^{y ~ ■),v\y\y') By^'^r^). 

yev y'&v\y 

It is worth pointing out that these hmits do not depend on ^. Hence, we have 
point-wise hmits for H^ron^ 



(Lv.aG)(r;) := -m\r,\G{r,) +Y. I Gi^U x)B]^-Hv\Od: 



X. (4.7) 



The convergences (|4.5p and (|4.6p in the space Cc are estabhshed by our next 
Lemma. 

Lemma 4.3. Let conditions of Lemma \4-l\ hold. Then, for a. a. x £ M."^ and 
for X-a.a. ^ G To, the convergence (|4.5p and (|4.6p hold in the sense of norm of 
Cc- 
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Proof. By dm and (03]), it is easy to see that e"!''! (i^o'^^e.B (2:,^ U •)) (v) has 
the form Ae(?7) + eB^irf). Moreover, the proof of Lemma [4.11 assures that i?£ 
has an integrable majorant. Hence, by the dominated convergence theorem, 
eB^ — !> in Cc- Next, using again the dominated convergence theorem and 
taking into account (|4.5|) and (|4.6p . we wiU be able to show convergence of 
to in Cc once we find uniform in e integrable estimate for the corresponding 
differences |Ae 

Since e"^"^ < 1 and ipeix) < 4'{x), for the establishment case, we have 
K+ ^ a+(x - y') 

y'erj 

g-eiJ*(x,c)g-s0(x-,')e^ (^^(^ _ ^ \ y') _ (-0(.T -•),'? \ 2/') 
+ ^ ^ a+(a;-2/)fe+(y-2/') 

-ex {-4>{x - O-'yV {2/,y'}) 

E a+(a;-2/)?^+(2/-2/')eA(0(x--),'7\{2/,y'})- 

ye-n y'en\y 

The last expression is an element of Cc, in view of (j2.9p and (|2.8p . Indeed, 



and, a similar equality holds for the second term. 

One can get the same result for the fecundity case in a similar way. □ 

Let us denote by the closed ball of radius c > in the Banach space 
L°°{R'^). 

Using Lemma [4.31 one can easily pass to the limit in (|4.ip . Therefore, in 
view of the general results presented in [M] we are able to state now the main 
theorem of this section. 

Theorem 4.4 (Proposition 4.2, Theorem 4.4. of [2]). Let the conditions of 
Lemma\4-1\ hold. Then 



1. (iy.u,^) are generators of a holomorphic semigroups Uv,i{t) on Cc. 

2. U^,f^(t) — > Uv,(i{t) strongly in Cc uniformly on finite time intervals. 
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3. There exists £ (0; 1) such that for any a G (ao; 1) the operator Ly"'^ 
Ly^ with the domain 

Dom(L©'J) = {A: e I ty^k G :^}. 

will he a generator of a strongly continuous semigroup Uy'^ (t) on the space 
K-ac ■ Moreover, for k G fCaC 

iL*y^k){7j) = -m\r^\k{7^) + J2 [ fc(eu(77\x))i?J'«(OdA(e). 
4- Let a G (ao; 1), po G -B^- Then the evolution equation 



it=o 



has a unique solution kt — ex{pt) in K-aC provided pt belongs to and 
satisfies the Vlasov-type equation 

^ft (x) = -mpt{x) + ^ exipu OB^'HOdX(.0- (4-8) 

Taking into account the explicit expressions for B^f^", one can rewrite (|4.8p 
in more simple form. Namely, using (|2.9p , for the establishment case we obtain 

^ Pt{x) = -mpt{x) + [ [ ex{pt,V^y)^~^ex{-(j){x - ■),7j)a+{x - y)dyd\{r]) 



dt 

" u " i^ -^ 

ex{pt,V U {y, y'})a+{x - y)h+{y ~ y')ex {-(j){x - •), v) dydy'dX{r]) 



and, by (|2.8p . we will have 
d 

g^Pt{x) = -mpt{x) + >c+{pt*a'^){x)cxp{-{pt*(t>){x)} (4.9) 

+ {{iPt * b+)pt} * a+)(x) exp{-(pt * (t>)ix)} . 

Here and below * means usual convolutions of functions on M.'^. 
Analogously, for the fecundity case, we obtain 

d 

-Q^Pt{x) = ~mpt{x) + K+{{ptexp{-pt* (j))} * a+){x) (4.10) 

+ {{{Pt *b+)pt exp{-pt *(f>)} *a+){x). 

Of course, we are mostly interesting in nonnegative solution of Vlasov equa- 
tion to have kt = ex{pt) is a correlation function of Poisson non- homogeneous 
measure with intensity pt- The existence and uniqueness of such solution we 
establishes by the following propositions. 
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Proposition 4.5. Suppose there exists A > such that < niax{a+(a;), b'^{x)} < 

ana 

x+(l + ^(<?i)) +c(6+) (2 + ^(0)) < m, (4.11) 

+ < m. (4.12) 

e 

Then the equation (|4.9I) mf/i initial < po & has a non-negative solution 
pt- Moreover, pt G B"^ and it is a unique solution from B^ . 

Proof. Let us fix some T > and consider the Banach space Xt = C([0; T],L°°{R'^)) 
of all continuous functions on [0; T] with values in L°°(M'^); the norm on Xt is 
given by 

IIuIIt := max \\ut\\ l-- (R'I) ■ 

te[0;T] ^ ' 

We denote by X^ the cone of all nonnegative functions from Xt- Denote also 
by B^ ^ the set of all functions u from X^ with ||w||t < c. 

Let $ be a mapping which assign to any v G Xt the solution ut of the linear 
Cauchy problem 

-g^ut{x) = -mut{x) + x+{vt * a+)(a;)exp{-(ut * (t)){x)] 
' + ({(wt * h+)vt} * a+)(x)exp{-(ut * 4>){x)], 

for a. a. a; G M''. Therefore, 

($v)t(a;) = e-"*po(2^) (4.13) 

+ /" e-™(*-'''>f+(u, *a+)(a;)exp{-(w, *,/))(2;)}ds 
Jo 

+ / e-™(*-"-)({(w«*6+)z;J*a+)(a;)exp{-(w«*(/.)(x)}ds. 
Jo 

It is easy to see that $w G X^. Indeed, one can estimate 

\{'^v)t{x)\ < \po{x)\ + {x+\\v\\T + {b+)\\vfT) f e-^'''^"'ds 

^ , ^+\\v\\T + {b+)\\vrT 
— c + , 
m 

where we have used the trivial inequality 

II/*5||l==(k-) < ll/llLnK'')ll3lli~(K-)' / e Li(R^), 5 e i°°(M'^). (4.14) 
Clearly, ut solves f|4.9p if and only if u is a fixed point of the mapping $ : Xt — )■ 
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We have that v E impUes $w G . Next, for 
|($u)f(x) - {^w)t{x)\ 



< K 



/' 

^0 



{vs * a+)(x) exp{-(ws * (/>)(x)} 



ds 



{{{vs *b+)vs} *a+){x) exp{-(ti5 * <P){x)} 
-{{{ws * b'^)ws} * a'^){x) exp{-(ws * 4>){x)} 



ds. 



Taking into account (|4.14l) and obvious inequalities e ^x < e ^ for x > 0, 

|e~° — e~''\ < \a — b\ for a,b > 0, and, moreover, 

Ipe^" - qe-''\ < e'^lp - q\ + qe-''\e-^''-''^ - 1\ < e-"|p - q\ + qe-''\a - b\, 
for any a, b,p,q > 0, we obtain 



-m(t — s) 



{vs * a+)(a;) exp{-(t;s * 4>)ix)} 



ds 



-{ws * a'^){x) exp{-(ws * 4'){x)} 
e~m(t^s) 1^(1^^ _ ^ ^+)(^) exp{-(i;, * q^){x)} 
^{ws * a^)(.x) exp{-(ws * (l)){x)]{\vs - Ws| * (l)){x)]ds 



< X \\v — w\\ 
and, similarly. 



\\v-w\\t(i + -{(^)) f e-"'^'-'^ds<\\v-w\\T — (l + -{^)); 



{{{vs * b+)vs} * a+)(a;) exp{-(us * (l)){x)} 



'{{{ws * b^)ws} * a+)(a;)exp{-(u;s * 4'){x)} 



ds 



< 



+ {{{ws * 6+)ws} * a+)(a;)exp{-(ti;s * - Ws\ * d^^.lS) 

Using the bound 

{{{ws * b+)ws} * a+)(a;)exp{-(ws * (l)){x)] 
< \\ws * 6+||Loo(Rd)(u;s * a+)(a;)exp{-(ws * 



22 



we may continue to estimate (|4.15l) as follows 



< 



ds 



-m(t~s) I 



(Vs * b'^)Ws — (Ws * b'^)ws ds 



For \\v\\t < c, I|u;||t < c one can estimate this expression by 

2c\\v - w\\T{b+) + c{b+)^\\v - w\\T{(t>)) e"^'^'~'Us. 

Therefore, for v,w £ , \\v\\t < c, ||w||t < c 

W^v - ^w\\t < — fl + -(0)) \\v - w\\t + ^^^(2 + -{(j))) \\v - w\\t. 
mVe/ m \ e / 

Moreover, if po G and v e B^^^ then, by (l4TT3t . 



C.4 



1 



me 



provided (|4TT2|) holds. 

As a result, by (|4.1ip . (|4.12p . $ is a contraction mapping on the closed set 
B^ ^. Taking, as usual, i;^"^ = <i>"w(°\ n > 1 for v^^^ G B^ ^ we obtain that 

C B+^ is a fundamental sequence in Xt which has, as a result, a unique 
limit point v € Xt- Since B^ ^ is a closed set we have that v £ B^ ^. Then, 
according to the classical Banach fixed point theorem, v will be a fixed point of 
4> on Xt and a unique fixed point on B^ 



T,c- 



□ 



The same considerations may be applied to the Vlasov equation (|4.10p . To 
combine these results with statement of Theorem 14.41 we need additionally that 
(jiHj) . (O^ hold with c = aC. 
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